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Abstract. Recent work in the literature has studied fourth-order elliptic operators on 
p manifolds with boundary. This paper proves that, in the case of the squared Laplace 
^ ■ operator, the boundary conditions which require that the eigenfunctions and their normal 
- ^ , derivative should vanish at the boundary lead to self-adjointness of the boundary- value 

XI 

' problem. On studying, for simplicity, the squared Laplace operator in one dimension, on a 
closed interval of the real line, alternative conditions which also ensure self-adjointness set 
to zero at the boundary the eigenfunctions and their second derivatives, or their first and 
third derivatives, or their second and third derivatives, or require periodicity, i.e. a linear 
relation among the values of the eigenfunctions at the ends of the interval. For the first 
four choices of boundary conditions, the resulting one-loop divergence is evaluated for a 
real scalar field on the portion of fiat Euclidean 4-space bounded by a 3-sphere, or by two 
concentric 3-spheres. 
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1. Introduction 



The current attempts to develop quantum field theories of fundamental interactions have 
led to the consideration of fourth-order or even higher-order differential operators on closed 
Riemannian manifolds [1-5], or on manifolds with boundary [6, 7]. The analysis of the 
transformation properties under conformal rescalings of the background metric g leads, 
in particular, to the consideration of conformally covariant operators P, which transform 
according to the law 



if g rescales as g^^ — e^^g, m being the dimension of the Riemannian manifold which is 
studied. One of the physical motivations for this analysis lies in the possibility to use the 
Green functions of such operators to build the effective action in curved space-times [5] . 

Another enlightening example is provided by the ghost sector of Euclidean Maxwell 
theory in vacuum in four dimensions. The corresponding field equations are well known 
to be invariant under conformal rescalings of g. On the other hand, the supplementary 
(or gauge) conditions usually considered in the literature are not invariant under con- 
formal rescalings of g. Even just in flat Euclidean 4-space, conformal invariance of the 
supplementary condition is only achieved on making the Eastwood-Singer choice [8] : 



where Ac is the electromagnetic potential (a connection 1-form in geometric language). 
The preservation of Eq. (1.2) under gauge transformations of A^. 




(1.1) 







(1.2) 



fAc = Ac + Wcf 



(1.3) 



is achieved provided that / obeys the fourth-order equation 



□V = o 



(1.4) 



where ^ is the box operator composed with itself: 



0^ = VaV«V6V^ 
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In the corresponding quantum theory via path integrals, one thus deals with two indepen- 
dent ghost fields (frequently referred to as the ghost and the anti-ghost), both ruled by 
I 1^, which is a fourth-order elliptic operator, and subject to the following boundary con- 
ditions (hereafter. Vat = A^"Va denotes the covariant derivative along the inward-pointing 
normal N°' to the boundary): 

[s]dM = (1.5) 



Vjve 



= 0. (1.6) 

dM 



Remarkably, since one now deals with a fourth-order elliptic operator, it is insufficient to 
impose just Dirichlet or Neumann (or Robin) boundary conditions. One needs instead 
both (1.5) and (1.6), which are obtained from the following requirements: 

(i) Gauge invariance of the boundary conditions on Af, [6, 7]. 

(ii) Conformal invariance of the whole set of boundary conditions [7] . 

(iii) Self-adjointness of the operator (see section 2). 

Although it remains extremely difficult to build a consistent quantization scheme via path- 
integral formalism for the full Maxwell field in the Eastwood-Singer gauge (the gauge-field 
operator on At, perturbations being, then, of sixth order [6, 7]), the investigation of the 
ghost sector remains of considerable interest in this case. There is in fact, on the one 
hand, the need to understand how to quantize a gauge theory in a way which preserves 
conformal invariance at all stages (as we just said), and on the other hand the attempt to 
extend the recent work on conformally covariant operators [1-5] to the more realistic case 
of manifolds with boundary. 

Our paper begins, therefore, with a detailed derivation of the boundary conditions 
which ensure self-adjointness of the operator For simplicity, the analysis is limited to 
one-dimensional problems, but the key properties are not affected by this sort of simplifica- 
tion. Section 3 proves the strong ellipticity and self-adjointness of the resulting boundary 
value problem. Section 4 studies the one-loop properties of a scalar field on a portion of 
flat Euclidean 4-space bounded by a 3-sphere, when the field is ruled by the Q ^ operator 
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and is subject to the boundary conditions (1.5) and (1.6). Section 5 extends the analysis 
of section 4 to the part of flat Euchdean 4-space bounded by two concentric 3-spheres. 
Results and open problems are discussed in section 6, and relevant details are described in 
the appendix. 

,4 

2. Self-adjointness of the operator 

We are concerned with the squared Laplace operator acting on scalar fields on a flat 
Euclidean background, in the case when curvature effects result from the boundary only. 
Moreover, motivated by quantum cosmology and Euclidean quantum gravity, the boundary 
is assumed to be a 3-sphere of radius a, or a pair of concentric 3-spheres [6, 7]. The former 
case, in particular, may be viewed as the limiting case when the wave function of the 
universe is studied at small 3-geometries (i.e. as a —> 0), as shown in [9]. 

In our problem it is hence possible to expand the scalar field on a family of 3-spheres 
centred on the origin, according to the familiar relation [10] 

oo 

e{x,T) = Y,en{T)Q^''\x) (2.1) 

n=l 

where r e [0, a], Q*-"^ are the scalar harmonics on a unit 3-sphere, S^, and x are local coor- 
dinates on S^. Thus, one is eventually led to study a one-dimensional differential operator 
of fourth order, and this makes it clear why all the essential information is obtained by 
the analysis of the operator B = on a closed interval of the real line, say [0, 1]. The 
operator B is required to act on functions which are at least of class (see (2.22)), and 
the following definition of scalar product (anti-linear in the first argument) is considered: 

{u,v) = I u*{x)v{x)dx. (2.2) 
Jo 

We now want to study under which conditions the operator B is sclf-adjoint, which means 
that it should be symmetric, and its domain D{B) should coincide with the domain of 
the adjoint B^ . For this purpose, we first study the relation between the scalar products 
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(Bu, v) and (u, Bv). We have then to integrate repeatedly by parts, using the Leibniz rule 
to express 



dx \ dx^ J ^ dx \ dx'^ dx J ^ dx \ dx dx^ J ^ dx \ dx^ J 



This leads to 



{Bu, v) 



d^ 



dx^ 

^ d^v 



n 1 



1 



d^u* dv 
dx"^ dx 



n 1 



+ 



du* d'^v 
dx dx"^ 



u 



dx^ 



+ {u,Bv). 



(2.3) 



Bearing in mind that the adjoint, B\ of ^ is again the operator it is thus clear that 
the condition {Bu,v) = {u,B^v) is fulfilled provided that both u e D{B) and v e D{B^) 
obey the same boundary conditions, for which the four terms expressing the difference 
are found to vanish. Some ways to achieve this are as follows. 



(i) First option: 



w(0)=u(l) = •u'(O) = = 
v{0) = v(l) = v'(0) = v'(l) = 0. 



(2.4) 
(2.5) 



(ii) Second option: 



u{0) = u{l) = u"{0) = u"{l) = 
v(0) = v(l) = ^;"(0) = -^"(1) = 0. 



(2.6) 
(2.7) 



(iii) Third option: 



u'{0) = u'{l) = u"'{0) = u"'{l) = 
^;'(0) = ^;'(1) = ^;'"(0) = v"'{l) = 0. 



(2.8) 
(2.9) 



(iv) Fourth option: 



u"{0) = u"{l) = u"'{0) = u"'{l) = 



(2.10) 
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^;"(0) = ^;"(1) = v"'iO) = v"'{l) = 0. 



(2.11) 



(v) Periodic boundary conditions: 



Ml) 
uiO) 

u'{0) 

u"{0) 7* 

^^^^(1) _ _1_ 
u"'{0) ~ j3* 



= 7 



1 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



and the same for v E D{B^)^ where (3 and 7 are some constants, not necessarily equal. 
This is made possible by the fourth-order nature of our operator. By contrast, if we were 
studying the first-order operator on the set of absolutely continuous functions on [0, 1], 
periodic boundary conditions leading to self-adjointness would involve one and the same 
complex parameter [11]. 

The solutions of the eigenvalue equation for the operator S, i.e. 



Bu 



= Xu 



(2.16) 



read 



u{x) = Ci cos px + C2 sin px + C3 cosh px + C4 sinh px (2-17) 

where p = A^/^. In particular, the periodic boundary conditions (2.12)-(2.15) lead to a 
linear algebraic system for the evaluation of the coeflBcients Ci, C2, C3 and C4 which admits 
non-trivial solutions if and only if the determinant of the following matrix vanishes: 



^ (cos p-P) 
— sin p 

(^-cosp+^^ 



\ 



sm p 



sm p 
(cos p — 7) 
— sin p 

- cos p + ^ 



(cosh p — (3) 
sinh p 

(^coshp- 
sinh p 



sinh p \ 
(cosh p — 7) 

sinh p 
(coshp- ^) / 
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The above determinant, denoted by S, turns out to have the form 



6 = Fi + F2(cos p + cosh p) + F3(cos p)(cosh p) 



(2.18) 



where 




(2/?7 + l) 



+ 2/?7 



(2.19) 




(/3 + 7) 
/?*7* 



(2.20) 



^3 = 2 



(/j + 27) (2/j + 7) 
13* 7* 



/3*7* 



1 



(2.21) 



To sum up, if the conditions (2.4) and (2.5), or (2.6) and (2.7), or (2.8) and (2.9), or 
(2.10) and (2.11), or (2.12)-(2.15) are satisfied, the domains of B and of its adjoint turn 
out to coincide: 



With our notation, ^C^[0, 1] is the set of functions in L^[0, 1] whose weak derivatives up 
to third order are absolutely continuous in [0,1], which ensures that the weak derivatives, 
up to fourth order, are Lebesgue summable in [0,1], and that all u in the domain are of 
class on [0,1]. Of course, symmetry of B is also obtained with the boundary conditions 
just described. 

In other words, at least five sets of boundary conditions, (i) or (ii) or (iii) or (iv) 

,4 

or (v), can be chosen to ensure self-adjointness of the operator Hereafter, we first 

consider the option (i), since, as was stated in the introduction, it is the one which agrees 

with boundary conditions motivated by the request of gauge invariance and conformal 

invariance, if the scalar field is viewed as one of the two ghost fields of Euclidean Maxwell 

theory in the Eastwood-Singer gauge. We also stress again that nothing is lost on studying 

just the "prototype" operator The one-dimensional fourth-order operator may take a 

more complicated form in some set of local coordinates (see (4.1)), but is always reducible 
,4 

to the operator ^ on the real line (more precisely, a closed interval of in our problems) . 



D{B) = D{B^) ^ { 



u:ue ^C"^[0,1], (2.4) or (2.6) 



or (2.8) or (2.10) or (2.12) - (2.15) hold} . 



(2.22) 
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3. Strong ellipticity and self-adjointness of the boundary value problem 

A number of points discussed in the previous section need to be put on firmer ground. In 
particular, we are concerned with the issue of eUipticity of the boundary value problem. 
This is studied in terms of the leading symbol of our differential operator, which is a squared 
Laplacian on a Riemannian manifold with smooth boundary. It is indeed well known that 
the Fourier transform makes it possible to associate to a differential operator of order k a 
polynomial of degree k, called characteristic polynomial or symbol. The leading symbol, 
(Jl, picks out the highest order part of this polynomial. For a squared Laplacian, denoted 
by F, it reads 

aL{F;x,0 = = g'^g'^'U-^Uai- (3.1) 

With a standard notation, x are local coordinates on M, and are cotangent vectors: 
Cm ^ T*{M). The leading symbol of F is trivially elliptic in the interior of M, since the 
right-hand side of (3.1) is positive-definite, and one has 

dim V 



/ \ / . \ aim V 

det(aL{F;x,0-X) =(\^\^-X) ^0 (3.2) 



for all A e C — R-|-, where V is the vector bundle over M whose sections are the physical 
fields (f, acted upon by F : C'^{V,M) C°°(y,M). In the presence of a boundary, 
however, one needs a more careful definition of ellipticity. First, for a manifold M of 
dimension m, the m coordinates x are split into m — 1 local coordinates on dM, denoted 
by {x''}, and r, the geodesic distance to the boundary. Similarly, the m coordinates are 
split into m — 1 cotangent vectors Q e T*{dM), jointly with a real parameter u e T*{R). 
The ellipticity we are interested in requires now that (tl should be elliptic in the interior of 
M, as specified before, and that strong ellipticity should hold. This means that a unique 
solution exists of the eigenvalue equation for the leading symbol: 

aL{F; {x''},r = 0, {C,} ,a; ^ -idr) - x\^{r) = (3.3) 

subject to the boundary conditions and to a decay condition at infinity. Before defin- 
ing these concepts, note that, in (3.3), iuj is eventually replaced by the operator of first 
derivative with respect to the geodesic distance to the boundary 



A complete formulation of boundary conditions needs some abstraction. For this 
purpose, one has to consider two vector bundles, Wp and W'p, over the boundary of M, 
with a boundary operator Bp, relating their sections, i.e. 

Bp : C°^{WF,dM) C°°{Wp,dM). 

All the information about normal derivatives of the fields is not encoded in Si? but in the 
boundary data i^ri'^) ^ C°°{WF,dM). For example, with boundary conditions involving 
(f and its first normal derivative, one has 

and the boundary conditions read Bpi/jpif) = 0. The sections of the bundle W'p, which 
remained unspecified so far, are obtained by applying to the sections of Wp the operator 
whose main diagonal coincides with the main diagonal of Bp. More precisely, if the 
boundary conditions are mixed, on writing Bp = PpL for some projector Pp : Wp Wp 
and some operator 

L : C°°iWp, dM) C°°{Wp, dM) 

one has ipp e C°°{Wp, dM) realized as ipp = Ppx, for some x £ C°°{Wp, dM). However, 
when the boundary operator (3.5) is considered, the projector Pp is turned into Bp, and 
the strong ellipticity condition demands that a unique solution of Eq. (3.3) should exist, 
subject to the boundary condition 

ag{Bp){{x''} , {Q})^Pp{^) = yi^M G C^{W^, dM) (3.6) 

and to the asymptotic condition 

lim ip{r) = 0. (3.7) 

r— >oo 
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With a standard notation [7, 12], ag{BF) is the graded leading symbol of the boundary 



Similarly to the case of the differential operator acting on physical fields, one is here 
mapping the boundary operator into its counterpart via Fourier transform. In the case of 
mixed boundary conditions for operators of Laplace type [13], Bp has off-diagonal elements 
which are first-order tangential operators, whereas complementary projectors occur on the 
main diagonal. One then finds a more elaborated structure [13]: 



where T is an anti-self-adjoint matrix. 

The asymptotic condition (3.7) picks out the solutions of the eigenvalue equation (3.3) 
which satisfy (3.6) with arbitrary boundary data '4>'f{v) ^^id vanish at infinite geodesic 
distance to the boundary. When all the above conditions are satisfied e T*(5M), VA e 
C - R+,V(C,A) ^ (0,0) and V?/'^ e C^{W'p,dM), one says that the boundary value 
problem {F,Bp) for the squared Laplacian is strongly elliptic with respect to the cone 
C — R-I-. Following [14], we find the solution of (3.3), (3.6) and (3.7) in the form (xi and 
X2 being some constants) 



operator Bp in the local coordinates {a:'^} , {Cj}- When Bp takes the form (3.5), ag{BF) 
may be defined by 




(3.8) 




^{r) = Xie + X2e 



-p2r 



(3.9) 



where, on setting 




(3.10) 



we define 




(3.11) 




(3.12) 
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Comparison with the case of the Laplace operator shows that one obtains strong eUipticity 
provided that ±-\/A G C — R-|-, which yields [14] 



Ae (C-R+)nC = C-R+. (3.13) 
The boundary condition (3.6) leads to the equation Ax = ip, where 

A=( M. (3.14) 

V-pi -P2J 

This matrix is trivially invertible, and hence existence and uniqueness of the solution is 
guaranteed. On writing ip' = where, according to the rule described after (3.5), 

one has (for some constants 71 and 72) 

A = 11 + 12 (3.15) 

tp'i = -pill - P2I2 (3.16) 

one finds 

XI = 71 (3.17) 

X2 = 72- (3.18) 

As stressed in [12], the condition of strong eUipticity is essential to ensure the existence 
of the asymptotic expansions normally assumed in the theory of heat-kernel asymptotics. 
In other words, if one cannot prove strong eUipticity for a given choice of boundary con- 
ditions, the local asymptotics of the fibre trace of the heat-kernel diagonal, and the corre- 
sponding, global asymptotics (resulting from integration over M) do not contain just the 
terms whose occurrence is ensured by invariance theory and by the Weyl theorem on the 
invariants of the orthogonal group [12]. There are, instead, highly singular contributions 
to the heat-kernel diagonal, so that their integral over M does not exist in any sense. 

It is therefore reassuring to see that, precisely in the case more relevant for quantum 
field theory [6, 7], strong eUipticity holds, and hence the resulting one-loop theory is 
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well defined. The following sections are devoted to a detailed evaluation of such one- 
loop approximation, but now we should clarify the self-adjointness issue for the squared 
Laplacian. For this purpose, let us recall that, given a symmetric operator A, with adjoint 
A^, its self-adjointness property can be studied by evaluating the dimension of the space 
of solutions of the equation A^ip = i(p, jointly with the corresponding dimension for the 
equation (f = More precisely, one defines the deficiency sub-spaces [11] 

n+ = Ker(z - A^) (3.19) 

n- = KeT{i + A^) (3.20) 
and the associated deficiency indices [11] 

n+(A) = dim[7^+] (3.21) 

n_(A) = dim[n-]. (3.22) 
Two theorems are then very useful [11]: 

Theorem 3.1. Given a closed symmetric operator A with deficiency indices n-|- and n_, A 
is self-adjoint if and only if n+ = = n_. Moreover, A has self-adjoint extensions if and 
only if n+ = n_ , and a one-one correspondence exists between self-adjoint extensions of A 
and unitary maps from onto H- . 

Theorem 3.2. If A is a symmetric operator with domain D{A), and if a conjugation C exists 
which maps D{A) into D{A) and commutes with A. CA = AC, then nj^{A) = n_(A), 
and hence A has self-adjoint extensions. 

,4 

In the case of the operator B = studied in section 2, since complex conjugation 
commutes with B, we immediately know from theorem 3.2 that the deficiency indices 
of B are equal. The solutions of the equations = on L^(0,oo) (we shall later 
restrict to [0,1]) are weak solutions. However, by virtue of the elliptic regularity theorem, 
these solutions are infinitely differentiable and hence strong solutions [11]. Now the strong 
solutions of 

= M^) (3-23) 
12 



can be written in the form (f{x) = e°'^, with a a root of the equation a'^ = i. One then 
finds 

ai = e*^ = cos ^ + i sin ^ (3.24) 
8 8 

a2 = ie* = — sin — \-i cos — (3.25) 
8 8 

as = — e* "s = — cos — — i sin — (3.26) 
8 8 

•TT TT TT / ^\ 

a4 = —le 8 — sm — — I cos — . (3.27) 
8 8 

Thus, only the strong solutions f2{x) = e"^^ and fsix) = e"^^ are in L^(0, cxo), and 
n+(i?) = 2. Similarly, the strong solutions of the equation 

ip(x) = -i(f{x) (3.28) 



can be written in the form (p{x) = e^^, with /? a root of the equation = —i. One then 
finds 

(3i = e"^^ = cos ^ - i sin ^ (3.29) 
8 8 

P2 = ■ie"*^ = sin — + i cos — (3.30) 
8 8 

/33 = -e"*t = -cos^ +isin^ (3.31) 

8 8 

/54 = -ie~'^ = - sin - - i cos - (3.32) 

8 8 

which implies that only the strong solutions <P3{x) = e^^^ and ^^{x) = e^"*^ are in L^(0, cxo), 
and hence n-{B) = 2 = n+(S). This property suggests also a non-trivial link with the 
strong ellipticity analysis, where we have seen that the asymptotic condition (3.7) selects 
only two of the original four contributions to the solution (3.9). 

However, on L^[0, 1], all strong solutions resulting from (3.24)-(3.27) and (3.29)- 
(3.32) are acceptable, and the deficiency indices of the operator B are, therefore, n+(i?) = 
n-{B) = 4. The domains for B studied in section 2 and summarized in equation (2.22) 
are hence correctly interpreted as domains of self- adjoint extensions of B. 
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4. One-loop divergence on the Euclidean 4-ball 

The definition and evaluation of functional determinants remains a topic of crucial im- 
portance in quantum field theory. Here the task is even more interesting, because we are 
studying a fourth-order elliptic operator on a manifold with boundary. As shown in [6, 7], 
the resulting eigenvalue equation for the modes occurring in the expansion (2.1) turns out 
to be, on the Euclidean 4-ball, 

rf4 ^ 6 rf3 (2^2 - 5) 



dr^ T dr^ dr"^ 



(2n2 + l) d (n2-l)2 



+ 



T 



3 dr t4 



£n = An£n- (4.1) 



Thus, on setting M = A^^, the solution of equation (4.1) is expressed by a linear combi- 
nation of Bessel functions and modified Bessel functions [7] , i.e. 

T T 

^^3,„^ + ^,„TO). (4.2) 



Since the Euclidean 4-ball consists of a portion of flat Euclidean 4-space bounded by a 
3-sphere, the coeflBcients ^2,n and ^4,n have to vanish Vn > 1, to ensure regularity of 
at the origin. One is thus left with scalar modes of the form 

£n(T) = \- A^^n • (4.3) 

T T 

We focus on a ^(0) calculation for such a set of massless modes, subject to the boundary 
conditions (see (1.5), (1.6) and (2.4), (2.5)) 

[en]dM = (4.4) 

[dSn/dT]9M = (4.5) 
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because the resulting 1-loop analysis remains crucial in the course of studying quantum 
theory as a theory of small disturbances [15] of the underlying classical theory. Our cal- 
culation relies on the technique developed in [16] and applied several times by the present 
authors (see [17] and references therein). The starting point is the remark that, since 
(■-functions are L^-traces of complex powers of elliptic operators, they admit an integral 
representation with the help of the Cauchy formula. For example, for a function / analytic 
in the domain bounded by a curve 7, one has 



where F is a function having zeros at zi,...,zl with multiplicities ni,...,nL, respectively. 
Thus, on choosing f{z) = z~'^, one finds the desired integral representation of the C,- 
function in the form 



where 7 is the contour in the complex- 2; plane which encircles all roots of the equation 
Q{z) = 0, with Q the function expressing the equation obeyed by the eigenvalues by virtue 
of the boundary conditions. The contour 7 is then deformed into a new contour 7, which 
encircles the cut in the complex plane of the function , coinciding with the negative 
real axis. After some technical steps, one eventually finds 



where I{—z,s) is the regularized infinite sum defined by I{—z,s) = ^(^) logQ(2;). 
More precisely, on denoting now by fn the function occurring in the equation obeyed by 
the eigenvalues by virtue of the boundary conditions, after taking out false roots (e.g. 
a; = is a false root of the equation J,^{x) = 0), and writing d{n) for the degeneracy of the 
eigenvalues parametrized by the integer n, one defines the function 




z)— log F{z) dz 



(4.6) 




(4.7) 




(4.8) 



00 




(4.9) 



n—no 
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What is very useful is the analytic continuation "/(M^, s)" to the complex-s plane of the 
function I{M^, s), which is a meromorphic function with a simple pole at s = 0, i.e. 



"/(M2, sY = + /^(M^) + 0(s). (4.10) 

s 

The function /poie is the residue at s = 0, and makes it possible to obtain the ^(0) value 
as [16] 

C(0) = /log + /pole(M2 = OO) - /pole(M2 = 0) (4.11) 

where Iiog is the coefficient of the log(M) term in as M — > oo. The contributions Iiog 
and /poie(oo) are obtained from the uniform asymptotic expansions of basis functions as 
M ^ OO and their order n ^ oo, whilst /poie(O) is obtained by taking the M ^ limit of 

the eigenvalue condition, and then studying the asymptotics as n — oo. More precisely, 
/poie(oo) coincides with the coefficient of ^ in the asymptotic expansion as n ^ oo of 

^d{n) log[poo(^)] 

where Pooin) is the n-dependent term in the eigenvalue condition as M ^ oo and n ^ oo. 
The /poie(O) value is instead obtained as the coefficient of ^ in the asymptotic expansion 
as n ^ OO of 

id(n) log[po(n)] 

where pQ{n) is the n-dependent term in the eigenvalue condition as M — > and n — > oo [16, 
17] . Although such a technique was originally developed for second-order elliptic operators, 
it can be easily generalized to study our fourth-order operators, provided that one bears in 
mind that the eigenvalues have now dimension [length]"^ (see notation after (4.1)). Hence 
one should replace by in (4.9)-(4.11), but the final results remain unaffected. 

In our problem, the equations (4.3)-(4.5) lead to the eigenvalue condition (denoting 
by a the radius of the 3-sphere) 

detf -^^^^""^ Vo (4 12) 

V -In{Ma) + MaI'^{Ma) -Jn{Ma) + MaJ'^{Ma) ) ^ ' 
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which guarantees that non-trivial solutions exist for the coefficients Ai^n and A^^n in (4-3). 
At this stage, on using the limiting form of Bessel functions and Jn when the argument 
tends to zero (see (A.1)-(A.4)), one finds that the left-hand side of (4.12) is proportional 
to M^"^ as M — > 0. Hence one has to multiply by M~^" to get rid of false roots. Moreover, 
in the uniform asymptotic expansion of Bessel functions as M — > oo and n — > oo, both / 
and J functions contribute a factor (see (A. 5), (A. 7) and (A. 12)). These properties 
imply that /log takes the value 

oo 
n=l 

The calculation of /poie(oo) relies on the asymptotic expansions (A. 5), (A. 7) and (A. 12) of 
the appendix. One then finds that no n-dependent term occurs in the eigenvalue condition 
(4.12), which implies 

/pole(oo) = 0. (4.14) 

Last, /poie(O) is obtained after working out po{n) for (4.12). For this purpose, we remark 
that, as M — > and n —>■ oo, the first line of the matrix in (4.12) consists of two elements 
both equal to ^'{n+i) ' whereas the second line consists of two elements both equal to 
(bearing in mind that all powers of (Ma/2) can be safely omitted, if one is interested in 
Po{n)). Hence one finds exact cancellation of the contributions to /poie(O), i.e. 

/pole(0)-0. (4.15) 

By virtue of (4.11) and (4.13)-(4.15) one finds, on the Euclidean 4-ball, 

C(0) = -ji; (4.16) 

for a real, masslcss scalar field. 

It is interesting to notice that the effect of cancellation of contributions to /poie(O) 
arises because, whilst recurrence formulae for derivatives of Bessel functions: 

2J^{z) = Jn-liz) - Jn+liz) (4.17) 
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and those for modified Bessel functions: 



(4.18) 



have different signs before Jn+i{z) and In+i{z), respectively, only J^-i^z) and In-i{z) 
survive in the limit as M — > 0, and hence the determinant (4.12) is equal to zero in that 
limit. On differentiating the recurrence formulae (4.17) and (4.18) the appropriate num- 
ber of times, one can easily check that the corresponding determinants for the boundary 
conditions (ii), (iii) and (iv) of section 2 are also equal to zero at M = 0. Thus, the 
contributions to /poie(O) with these boundary conditions vanish as well. 

Moreover, the expression (4.13) for /log is trivially modified by adding to (— 2n) the 
integer numbers 1 for (ii), 3 for (iii) and 4 for (iv), respectively. Their contributions are 
proportional to 



Thus, the results (4.13) and (4.16) hold for the first four types of boundary conditions. 
5. One-loop divergence in the two-boundary problem 

In the two-boundary problem one studies a portion of flat Euclidean 4-space bounded by 
two concentric 3-spheres. This case is very interesting because it is more directly related 
to the familiar framework in quantum field theory, where one normally assigns boundary 
data on two three-surfaces (it should be stressed, however, that unlike scattering problems 
we are considering a path-integral representation of amplitudes in a finite region). 

On denoting by a and 6, with a > b, the radii of the two concentric 3-sphere boundaries, 
we can consider the complete form (4.2) of our scalar modes, because no singularity at the 
origin occurs in the two-boundary problem, and hence all linearly independent integrals 
are regular, for all r e [b,a\. We now impose the boundary conditions (1.5) and (1.6), 
which lead to the eigenvalue condition 



CO 



^n2 = Ci^(-2) = 0. 



n=l 



det 



/ IniMb) K^iMb) Jn{Mb) Nn{Mb) \ 

F,„(M6) FK^Mb) Fj^iMb) F^^Mb) 

In{Ma) Kn{Ma) MMa) iV^(Ma) 

\Fj^{Ma) FK^Ma) Fj^Ma) F^^Ma)) 



= 



(5.1) 
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where, for Z = I,K,J or N, we define 



Fz^iMx) = -Zn{Mx) + MxZ'^{Mx). (5.2) 

On using the approximate forms (A.1)-(A.4) of Bessel functions, one finds that the left- 
hand side of (5.1) is proportional to M° as M — > 0. Hence there are no false roots 
of (5.1). As a next step we notice that, as M ^ oo and n — >^ oo, the 1^ K^, Jn 
and Nn functions contribute a factor, whereas Fz„, defined in (5.2), contributes a 

y/M factor. Moreover, the dominant contribution to (5.1) as M — > oo and n — > oo 
is given by Kn{Mb), Nn{Mb), FK^{Mb), FN^{Mb) (in the first two rows), jointly with 

In{Ma), Jn{Ma), Fi^{Ma), Fj^{Ma) (in the last two rows). One then finds that I\og van- 
ishes, because 

1 °° 

n=l 

The value of /poie(oo) also vanishes, because no n-dependent term occurs in equation (5.1) 
when n — > oo and M — > oo. Moreover, /poie(O) vanishes as well, since the determinant 
leading to po{n) takes the form (by virtue of (A.1)-(A.4)) 



D{n) = det 



/ ir(n) -ir(n) \ 

-ir(n + l) ^r(n)(n + l) 

r(n+i) " r(n+i) ^ 

V r(n+l) ^ r(n+l) ^ / 



(5.4) 



and this vanishes exactly. 

To sum up, we find that the ({0) value is zero in the two-boundary problem: 

C(0) = 0. (5.5) 

With the same arguments presented in the end of section 4, the result (5.5) is found to 
hold for the first four boundary conditions studied in section 2. 



19 



6. Concluding remarks 

The analysis of the squared Laplace operator in flat Euclidean backgrounds is motivated 
by the ghost sector of Euclidean Maxwell theory in a conformally invariant gauge, but 
has been here restricted to a real scalar field. Further motivations result from the theory 
of conformally covariant operators, which is an important branch of spectral geometry, 
and finds applications also in Euclidean quantum gravity [1-7]. The contributions of our 
investigation are as follows. 

(i) The boundary conditions for which the squared Laplace operator is self-adjoint have 
been derived (cf [18]), taking as prototype the operator ^ on a closed interval of the real 
line. Interestingly, at least five sets of boundary conditions are then found to arise, and the 
option described by (2.4) and (2.5) coincides, if the field in (2.1) were a ghost field, with 
the boundary conditions obtained from the request of gauge invariance of the boundary 
conditions on A^, when the Eastwood-Singer supplementary condition is imposed. The 
general reader, however, should be aware that the above boundary conditions have already 
been studied in the mathematical literature. For example, Eq. (1.5.51) of [18] studies the 
more involved boundary value problem 

(a2+AA)v = / in Q (6.1) 

v = at dQ (6.2) 

_ 1 

Vnv = (l - Ar^ \ ai dVL (6.3) 

where Ar is chosen in such a way that ^1 — Ar^ ^ is a suitable bijective operator [18]. 

Here it is enough to remark that, since A2 + AA > A2 for all A G R_, the solution of 
(6.1)-(6.3) is unique for all smooth data. Hence strong ellipticity follows [18]. 

(ii) Given the fourth-order eigenvalue equation (4.1), the contribution of the corresponding 
eigenmodes to the one-loop divergence has been derived for the first time on the Euclidean 
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4-ball (see (4.16)), or on the portion of flat Euclidean 4-space bounded by two concentric 
3-spheres (see (5.5)). 

In our opinion, the property (i) is crucial because no complete prescription for the quan- 
tization is obtained unless suitable sets of boundary conditions are imposed, and sections 
2 and 3 represent a non-trivial step in this direction. The result (ii) is instead relevant 
for the analysis of one-loop semiclassical effects in quantum field theory. In other words, 
if one has to come to terms with higher order differential operators in the quantization of 
gauge theories and gravitation, it appears necessary to develop techniques for a system- 
atic investigation of one- loop ultraviolet divergences, as a first step towards a thorough 
understanding of their perturbative properties. 

Some outstanding problems are now in sight. First, it appears interesting to extend our 
mode-by-mode analysis to curved backgrounds with boundary. In this case, the fourth- 
order conformally covariant differential operator is more complicated than the squared 
Laplace operator, and involves also the Ricci curvature and the scalar curvature of the 
background. Second, one should use Weyl's theorem on the invariants of the orthogonal 
group to understand the general structure of heat-kernel asymptotics [12] for fourth-order 
differential operators on manifolds with boundary. A naturally occurring question within 
that framework is, to what extent functorial methods [7, 12] can then be used to compute 
all heat-kernel coefficients for a given form of the differential operator and of the boundary 
operator. Third, the recently considered effect of tangential derivatives in the boundary 
operator [13, 19-21] might give rise to generalized boundary conditions for conformally 
covariant operators. The appropriate mathematical theory is still lacking in the literature, 
but would be of much help for the current attempts to understand the formulation of 
quantum field theories on manifolds with boundary [7, 13, 17, 19-21]. 
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Appendix 

In sections 4 and 5 we need the asymptotic expansions of Bessel functions when the argu- 
ment tends to zero, or when both the argument and the order are very large [22]. In the 
former case, one finds, for all n > 1, 

(A.1) 

Kn(x) ~ ir(n)(a^/2)-» (A.3) 
JV„W~-ir(n)(x/2)-». (A.4) 

TT 

Moreover, when the argument is greater than the order, both being large, one finds, for 
modified Bessel functions [23] , 

/.M^ (A.5) 



Kn{nz) ~ E(-l)^^ (A.6) 



* s=0 



s=0 



where 



^v/m^ + log— -4==. (A.9) 
1 + Vl + -z'' 
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We do not need, in our calculations, the detailed form of the polynomials Ug and Vg, which 
are generated by recurrence relations and can be found in [23]. For the Bessel functions 

Jn{nz) and Nn{nz), with both n and z very large, we need the asymptotic expansions in 
section 8.41 of [22], which are more conveniently expressed after setting z = sec/3: 



Jn{n sec P) 



mr tan P 



cos tan /3 — n/3 — ^ 

m=0 



+ sin (n tan 3 — nd — ^ — -^-^ ^"^"'"o^.i 



(-l)"^r(2m + ^) A2m 

r(|) (ftan/?)^^ 

(A.IO) 



m=0 



Nn{n sec (3) 



nn tan /3 
cos tan /? — n/3 — 



sin n tan 



(-irr(2m+i) A2„. 



- m=0 ^(2) 

(-l)-r(2m+|) ^2^+1 



m=0 



rQ) (ftan/3) 



2m+l 



(f tan/3)'"^ 
(A.ll) 



where and A2m+i are numerical coefficients. In particular, it is crucial to pick out 
the dominant terms of the expansions (A.IO) and (A.ll), i.e. 



Jn{nz) ^ \ —{z^ - 1) 4 
mr 



cos n ^J z^ — 1 — arccos - 



TT 



+ 0(n-^) 



(A.12) 



Nninz) 



-iz^-iy 

mr 



sm n 



arccos - 
z 



+ 0(n-i) 



(A.13) 



Thus, in (4.12), both Jn{Ma) and J'niMa) contribute at large M (as weU as In{Ma) 
and I^{Ma)). In section 5, all Bessel functions and their first derivatives contribute, for 
the same reason, a factor ^Jj^ in the eigenvalue condition when both n and M tend to 00. 
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